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We propose a model for the liquid- liquid (Lq, L^f ) phase transition observed in osmotic pressure 
measurements of certain charged lamellae- forming amphiphiles. The model free energy combines 
mean-field electrostatic and phenomenological non-electrostatic interactions, while the number of 
dissociated counterions is treated as a variable degree of freedom that is determined self-consistently 
The model, therefore, joins two well-known theories: the Poisson-Boltzmann theory for ionic solu- 
tions between charged lamellae, and Langmuir-Frumkin-Davies adsorption isotherm modified to 
account for charged adsorbing species. Minimizing the appropriate free energy for each inter lamel- 
lar spacing, we find the ionic density profiles and the resulting osmotic pressure. While in the 
simple Poisson-Boltzmann theory the osmotic pressure isotherms are always smooth, we observe 
a discontinuous liquid-liquid phase transition when Poisson-Boltzmann theory is self-consistently 
augmented by Langmuir-Frumkin-Davies adsorption. This phase transition depends on the area 
per amphiphilic headgroup, as well as on non-electrostatic interactions of the counterions with the 
lamellae, and interactions between counterion-bound and counterion- dissociated surfactants. Cou- 
pling lateral phase transition in the bilayer plane with electrostatic interactions in the bulk, our 
results offer a qualitative explanation for the existence of the Lq, L^' phase-transition of DDABr 
(didodecyldimethylammonium bromide), but its apparent absence for the chloride and the iodide 
homologues. More quantitative comparisons with experiment require better understanding of the 
microscopic basis of the phenomenological model parameters. 



I. INTRODUCTION 

From naturally occurring phospholipids to synthetic 
double-chain surfactants, over a wide range of concen- 
trations, amphiphiles in aqueous solutions self- assemble 
into multilamellar phases. The stability of the lamellar 
stack depends not only on the type of amphiphile, but 
also on the competition between several inter-lamellar in- 
teractions Q. Attractive van der Waals interactions are 
balanced by repulsive interactions. Hydration repulsion 
usually dominates when the intervening water layer spac- 
ings are small (typically < Inm) or intermediate, while 
electrostatic and 'steric' undulation interactions usually 
prevail at intermediate to large spacing, up to hundreds 
of nanometers 0, 0, 0, IS • 

Between charged surfactants, the stabilizing repulsion 
is typically provided by the strong Coulomb interaction 
mediated by dissolved counterions and salt 0. These 
interactions are particularly strong for lamellar-forming 
charged surfactants whose counterions fully dissociate 
into solution. Salt can attenuate such electrostatic in- 
teractions via ionic screening. For surfactants that form 
flexible layers, the complicated yet important coupling 
between layer elasticity, undulations, and electrostatic in- 



teractions must also be considered 0, H, IS • But even 
when the effects of layer flexibility can be ignored, elec- 
trostatic interactions in multilamellar charged systems 
are non-trivial and, in general, difficult to understand 
because of the intimate link between counterion dissoci- 
ation, ionic screening, and ion-specific non-electrostatic 
interactions [TlLIT^ Il^. 

The determination of the collapse pressure of mem- 
brane stacks by Dubois et al. 14] adds a new 
twist. Synthetic cationic double-chain surfactant, dido- 
decyldimethylammonium (DDA+) with bromide as coun- 
terion (DDABr), is used to form a thermodynamically 
stable lamellar phase. The system undergoes a phase 
transition from a swollen liquid-like (Lq,) lamellar phase 
to another, more condensed, liquid-like lamellar phase 
(Lq/). This phase transition is induced by externally 
applied osmotic pressure; it is seen as a plateau in the 
osmotic pressure versus inter-lamellar spacing isotherms. 
Measured by small angle X-ray scattering, the abrupt 
change in spacing is typically between lOA to lOOA. In 
contrast, for the same surfactant with the bromide coun- 
terion replaced by chloride, DDACl, there is no evidence 
of a first-order transition [I^- fact, the experimental 
isotherm can be well fit by the usual Poisson-Boltzmann 
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(PB) theory 0. Further, with an iodide counterion the 
stack made of DDAI surfactant remains collapsed and 
does not swell at all 15]. Remarkably, a discontinuous 
increase in area per surfactant with no discernible in- 
plane positional order was experimentally found to co- 
incide with the collapse in bilayer spacing. Clearly, this 
collapse is strongly coupled to a lateral rearrangement 
in the bilayer plane, and cannot be solely the result of 
neutralizing surfactant headgroups by their counterions. 

It should perhaps come as no surprise that different 
halide counterions interact differently with the charged 
DDA+ surfactant layers. The number of electrons, hence 
properties like polarizability, vary widely for these ions, 
and we expect that ion-membrane interactions will be 
different too. By ranking ions according to their effi- 
ciency in salting-out proteins from solution, Hofmeister 
was first to observe — over a century ago — that differ- 
ent ions partition differently at aqueous interfaces [T^ . 
The Hofmeister rankin g is sur p risingly insensitive to the 
details of the interface jT3,[lMll3- Often, however, the 
preferential interaction follows the size and polarizability 
of the ion; large ions tend to be less repelled from (or 
more attracted to) oily interfaces. 

It has been proposed that the added van der Waals 
attraction of the ions to the higher index of refraction 
material may explain the Hofmeister ranking [23, 0, 113, 
I23L [2^ . Water ordering around ions at the interface that 
is structured differently from the bulk can also discrim- 
inate between ions. More polar izable ions, for example, 
may be attracted to ordered water molecules at the in- 
terface because of favorable interaction between dipoles 
and induced-dipoles. Due to their amphiphilic, liquid- 
like nature, surfactants present a special and complex 
interface to water and salt ions. However, using mea- 
surements such as electrophoretic mobility, nuclear mag- 
netic resonance (NMR), and buoyancy density- matching, 
a Hofmeister-like rankin g of anions has also emerged for 
ions at lipid interfaces |25l [2^ . Also, for single 
chain micelle- forming cationic surfactants, the area per 
molecule follows the Hofmeister series, increasing more 
in the presence of larger ions. j2^ . 

As is evident from NMR experiments, different ions 
not only associate differently with the amphiphile-water 
interface, but their binding may also restructure the in- 
terface itself ^26]. Computer simulations indicate that 
the restructuring of the amphiphilic headgroup region 
should be strongly influenced by the counter- ion size [2^ . 
Such conformational changes at the interface are possi- 
ble sources of non-ideal lipid mixing, because ion binding 
at the interface may effectively create two incompatible 
types of lipids: ion-bound and ion-detached. For exam- 

1)le, both experiments and simulations of lipid mixtures 
30, 31, 32, 33, 34, 35] show that charged and uncharged 
lipids tend to demix so as to minimize the line tension be- 
tween the different mismatched lipid species, often lead- 
ing to lipid lateral phase separation in the membrane 
plane. Perhaps most compelling are the phase transi- 
tions from lamellar to inverted hexagonal phases of pure 



DOPS induced by varying pH that changes the fraction 
of charged to uncharged ionizable lipids 36]. 

Can an added non-electrostatic attraction of ions to 
the lipid-water interface explain the observed transition 
for Br~ ions? The charge regulation model of Ninham 
and Parsegian indicates that while an added attrac- 
tion can significantly modify pressure isotherms, it can- 
not account for a first-order phase transition. 

Here, we propose a phenomenological model that 
explains the first-order phase transition in terms of 
an added coupling between electrostatics and non- 
electrostatic specific interactions at the interface. The 
model is motivated by the experimentally observed 
lamellar-lamellar phase transition in charged surfactant 
systems |Q and is a relatively simple extension of 
Poisson-Boltzmann theory. 

The gist of our model is to consider the possibility that 
a fraction of the counterions are not dissociated from the 
lamellar-forming cationic DDA+ surfactant, but rather 
stay associated with it on the membrane plane to form 
a neutral complex. The degree of dissociation is taken 
as a variational parameter in our free-energy formula- 
tion, and is optimized for each inter-lamellar distance 
[T^ . Further, we consider each lamella as composed 
of a binary mixture of neutral (associated counterions) 
and charged (dissociated counterions) surfactant species. 
Assuming an effective attractive 2nd order virial coeffi- 
cient between the two species, we find possible lateral 
phase separation in the lamellar plane forming neutral- 
surfactant rich and charged-surfactant rich phases, much 
as in regular solution theory 0, [s^- Like any phe- 
nomenological model, our model relies on several param- 
eters whose exact molecular origin is not well known at 
present. However, using reasonable values of these pa- 
rameters, we are able to fit well the experimental data. 

Our model couples the Poisson-Boltzmann theory for 
the counterions in solution with the Langmuir-Frumkin- 
Davies adsorption model that regulates the amount of 
dissociated counterions 1^] • It is, there- 

fore, an extension of the charge regulation model of Nin- 
ham and Parsegian 0. Analogous coupling between 
surface transitions and bulk interactions has been ana- 
lyzed in the context of hydration forces as well 
as electrostatic interactions j44] . 

In the model, the differences between monovalent ions 
(Cl~, Br~, and I~) are accounted for by using differ- 
ent interaction parameters between ions and lipid inter- 
face and between ion-bound and ion-detached lipids. To- 
gether with the repulsive hydration force, known to act 
strongly at small separations such as those found in the 
collapsed phase, these interaction parameters are suffi- 
cient to reproduce the experimental observations. Our 
main emergent result supports a lamellar-lamellar phase- 
transition as function of externally applied osmotic pres- 
sure for ions such as Br~. At small osmotic pressure 
and large inter-lamellar distances, most of the Br~ ions 
are dissociated, and the isotherm follows the PB re- 
sult 05 El- However, for larger pressures and smaller 
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separations, a large fraction of the Br~ ions remain asso- 
ciated, causing a lateral phase transition. Because of the 
coupling between electrostatics and the entropy of ions 
in solution, the lateral phase transition also leads to the 
discontinuous jump in inter-lamellar spacing witnessed in 
the osmotic pressure isotherm. 

We further consider the effect of added salt on the 
equilibrium state of the system. The main effect of salt is 
to screen electrostatic interactions, to reduce the coupling 
between the layers, and thus to diminish the magnitude 
of the first-order transition jump Indeed, our model 
predicts that for more than a critical amount of salt the 
phase transition disappears altogether. 

The outline of the paper is as follows: in Sec. II 
we present an extension to the usual PB theory, taking 
into account the non-electrostatic degrees of freedom and 
treating separately the counterion-only and added-salt 
cases. In Sec. Ill we present our numerically calculated 
isotherms to show the possibility of a lamellar-lamellar 
phase transition. We then discuss the link with experi- 
ments and comment on ion- specific effects. In Sec. IV 
we discuss our findings and end in Sec. V with a short 
summary and remarks on possible future directions. 



II. EXTENDED POISSON-BOLTZMANN 
THEORY 

A. Model 

The lamellar stack is composed of bilayers of double- 
chain surfactants such as DBA separated by regions of 
aqueous solution. The charged surfactant hydrophilic 
headgroups point towards the water region, while the hy- 
drophobic tails are packed in the inner lamellar region, 
away from the polar water environment; for DDABr, the 
thickness of the hydrocarbon part of the bilayer is of the 
order of 24-26A. The lamellar stack can be modeled as 
a one-dimensional periodic system. This approximates 
the lamellar lateral extent as infinite and each lamella as 
perfectly planar and rigid. We consider only the unit cell 
of the lamellar stack, as is depicted in Figure 1. 

For convenience, the unit cell width D spans only the 
aqueous inter- lamellar region, while the periodicity of 
the lamellar stack includes also the bilayer thickness: 
D + Dm. All local quantities depend only on the per- 
pendicular coordinate, z. Because the solubility of sin- 
gle surfactants (like DBA) in water is extremely low, we 
assume that all surfactants reside within the lamellar bi- 
layers. These bilayers are the two bounding interfaces of 
the unit cell. 

The model system, hence, is composed of two charged 
interfaces separated by a distance D. The system is over- 
all electroneutral; the amount of counterions in the aque- 
ous region is exactly balanced by the amount of charged 
surfactant on the two interfaces. Note that because of the 
insolubility assumption, the ratio between the lamellar 
bilayer thickness and D uniquely determines the relative 




FIG. 1: A schematic view of the lamellar phase. The unit cell width 
includes only the water region, —D/2 < z < D/2, and is composed of 
two charged surfaces located at z = ±D/2 (each being one leaflet of 
the surfactant bilayer). The counterions can adsorb on the two surfaces 
or dissociate into the inter-membrane water region of thickness D. 



concentration of surfactant and water. 

The midplane of the unit cell, see Figure 1, is chosen 
at z = 0, from which the two interfaces bearing equal 
charge densities are located at z = -^D/2. Because the 
system is symmetric about the z = midplane, it is 
enough to consider a unit cell in the range < 2: < 
D/2. The counterion concentration (mole per liter) c{z) 
and the mean- field electrostatic potential il){z) depend on 
the perpendicular coordinate z. All variables calculated 
at the lamellar surface, z = D/2^ will be denoted by 
a subscript s (e.g., V^g), while those calculated at the 
symmetric midplane, z = 0, by a subscript m (e.g., '0m)- 

Our aim is to calculate the equation of state for this 
lamellar symmetry and to express it as a relation between 
the thermodynamic variables: inter-lamellar distance, D., 
and osmotic pressure, 11. Just as in the van der Waals 
phenomenological theory of phase transitions, we search 
for thermodynamically stable states or stable branches 
of the free energy. When we find more than one branch, 
we can use the Maxwell construction to obtain the coex- 
istence region. 

Minimization of the model free energy with respect 
to electrostatic and non-electrostatic degrees of freedom 
(see below) will eventually lead to the equation of state, 
our final goal. The overall free energy + + Fhyd 
per unit cell area is a sum of volume contributions com- 
ing from the electrolyte solution within the cell, including 
ion (electrostatic) Fy terms, surface contributions having 
their origin at the interfaces Fg, and hydration interac- 
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tion Fhyd, dominant at small D separations. 

Because we approximate the hydration interaction as a 
separable term, we can independently minimize the elec- 
trostatic contribution to the free energy. As we show 
next, these interactions already suffice to account for 
a discontinues phase transition, while Fhyd is needed 
only to account for the experimentally found pressure at 
small D. We therefore discuss only the minimization of 
^tot = -\- Fs below, and return to present and discuss 
the added Fhyd when we compare theory and experiment 
in the Results section. 

Because the system is extensive in the interfacial area, 
the ion (electrostatic) volume free energy per unit area, 
Fy, is taken as the appropriate intensive quantity 



Fy 



nD/2 

Jo 



-ec^- — {^y^kBTc{\n{c/co)-l) dz . 

(1) 

The first two terms include the standard electrostatic en- 
ergy, with £ :^ 80 the dielectric constant of water. The 
last term is the entropy of mixing in the dilute limit, 
where cq is a reference concentration, T the temperature, 
and kB the Boltzmann constant. Because all counterions 
in solution originate from a surfactant molecule, their 
integrated concentration (per unit area) must be equal 
in magnitude and opposite in sign to the surface charge 
density 



nD/2 
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c{z) dz 



(2) 



This is also a charge density condition and can be trans- 
lated via Gauss' law into the electrostatic boundary con- 
dition (in gaussian units): il)'{D/2) = ^jj'^ = Aiicj/e^ link- 
ing the surface electric field ?/;^ with the surface charge 
density cr. Unlike the usual PB theory where either the 
surface charge or the surface potential is held fixed, here 
cr is a self-adjusting parameter which will be determined 
variationally from minimizing the total free energy Ftot- 
The second part of the total free energy comes from 
the surface free energy contributions of the amphiphiles 
residing on the planar bilayers. The surface free energy 
Fg has electrostatic and non-electrostatic energy terms 
as well as a lateral mixing entropy contribution. Ex- 
pressed in terms of the surface area fraction rjs = o?(j/e 
of charged surfactants. 



a Fs = etjjsVs - OLr]s 



:XVs 



+ ksT [r]s \nr]s + (1 - r]s) ln(l - r]s)] , (3) 

where the first term couples between the surface charge 
and surface potential. The other terms are the enthalpy 
and entropy of a two-component liquid mixture: charged 
surfactant with area fraction 77^ and neutralized, ion- 
bound surfactants with area fraction 1 — rjs- The pa- 
rameters a and x phenomenological, respectively de- 
scribing the counterion-surfactant and the surfactant- 
surfactant interactions at the surface. As in the charge 



regulation model ^^^^ a < means that there is 
an added non-electrostatic attraction (favorable adsorp- 
tion free energy) between counterions and the surface; 
the more counterions are associated at the surface, the 
smaller the amount of remaining charged surfactant. 

The parameter x is the most crucial and unique ele- 
ment in our model. Representing non-ideal mixing ten- 
dencies in the bilayer plane as in regular solution theory, 
it alone (together with the usual components of stan- 
dard PB theory) is sufficient to account for a coupled 
transition in the bilayer plane and in the bulk. As in the 
Frumkin adsorption model js^, 113, 113, a posi- 

tive X parameter represents the tendency of surfactants 
on the surface to phase separate into domains of neutral 
and charged surfactants. 

Changing to dimensionless variables, we define y{z) = 
e^{z)/kBT, (/)(z) = a^c{z), a = a/ksT, x = x/ksT, 
and take for convenience coa^ = 1. Then, the total free 
energy is written as a functional of the variables y{z), 
(j){z), and a function of 7^5, and includes the conservation 
condition, eqI21 via a Lagrange multiplier, /i: 
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(4) 



Next, we minimize Ftot with respect to the surface 
variable 775, and the two continuous fields (j){z)^ vi^)' 
dFtot/drjs = SFtot/S(t){z) = 5Ftot/6y{z) = 0, correspond- 
ing to three coupled equations of state 
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(5) 
(6) 

(7) 



The free energy Ftot is also a function of surface potential 
ys and of the inter-lamellar spacing D. The differentia- 
tion with respect to D gives the osmotic pressure (to be 
discussed in Sec. II. B), while the variation with respect 
to ys gives the usual electrostatic boundary condition: 



y'{D/2) = y', 



(8) 



where Ib = e'^/{ekBT) is the Bjerrum length, equal to 
about 7A at room temperature for aqueous solutions 
{e = 80). The Lagrange multiplier, /i, acts as a chem- 
ical potential, but with the important difference that it 
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is related not to the bulk reservoir concentration, but 
rather to the concentration at the midplane, For 
a single counter ion type, we can choose, without loss 
of generality, the potential at the mid-plane to be zero, 
ym = et/jm/ksT = and then from eos 16171 



y"{z) 



47rl 



V(z) 



(9) 



Not surprisingly, we recover the Poisson-Boltzmann 
equation [31 connecting the electrostatic potential y{z) 
with the counterion concentration (j){z) in the solution. 
This can be expected since the non-electrostatic contri- 
butions enter only via the surface interactions expressed 
in eq[5J 

Rewriting eqISJ we arrive at an expression similar to 
the Langmuir-Frumkin-Davies adsorption isotherm 39] , 
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Q-a-xris+Vs 



(10) 



with the following modifications: the concentration (j)^^ 
at the midplane replaces the bulk concentration (because 
of the constraint of overall charge neutrality), and the 
surface interaction energy in the exponent contains the 
electrostatic part ys . The simpler Langmuir isotherm 
is recovered in the limit of non-charged surfaces, = 0, 
and no surface interaction, x = 0- 
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(11) 



Here, there is a unique relation between (prn and 775, while 
in the general case of non-zero ys and x-> generalized 
Langmuir-Frumkin-Davies equation offers a transcenden- 
tal relation between rjg-, 4>m and ys with the possibility of 
more than one solution. 

The solution of the PB equation |H1 for two symmetric 
charged surfaces separated by a distance each having 
a surface charge density of a, has a well-known analytic 
form |3,[iln 



y{z) = — In [cos^(i^z) 



(12) 



where the constant K is determined from the electro- 
static boundary condition y'{D/2) = y'^ = AttIbtIs/o?^ 
eqlHl as 

KDt8in(KD/2) = '^^rjsD . (13) 

We can now express ys and (j)m as function of a single 
dimensionless variable u = KD/2: 
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(14) 



with the constants Ci and C2 obtained as 



Ci = 


a2 


-kIb 


C2 = 








ttIb 



(15) 



The solution of the above equations, together with the 
adsorption isotherm eq ^1 completely determines the 
counterion density profile and the mean electrostatic po- 
tential via the solution K = K{D^ a, x)- 



B. Equation of State and the Lq 
Transition 



Phase 



We now solve the basic set of equations derived in 
the previous section. Substituting (pra-, ^s, and ys into 
the Langmuir-Frumkin-Davies isotherm, eq El we get 
an equation for u 



D 



Ciii tan(ii) 



1+ 



Co 



cos^(ii) 



f- 



j exp I — a — x-^wtan(w) 
(16) 

Obviously, for each imposed distance we can extract 
u = KD/2 from eq El Once u is known, it can be 
substituted into eq El wherefrom r^g, ys follow. 

Their values completely determine the potential profile 
iIj{z) = kBTy{z)/e and counterion profile c(z) = (j){z)/a^. 

If there is only one solution for the system has one 
stable state. Multiple u solutions indicate the possibility 
of coexistence between several phases, as well as transi- 
tions between them. Each phase corresponds to a sepa- 
rate branch of the free energy with its own dependence 
on D. If all solutions of eqElare non-zero, then we can 
have a first-order transition between two stable phases. 

The free energy as a function of D and r]s{D) can be 
obtained by substituting the results of minimization back 
into eqSl yielding 
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(17) 



The appropriate isotherm is now obtained by taking the 
derivative of the free energy with respect to giving 
the osmotic pressure Il{D) as 



n(Z)) = 



dFt. 



dD 



(18) 



As is usual in PB theory jj, the osmotic pressure can 
also be calculated from the contact theorem^ which relates 
the osmotic pressure with the value of the counterion 
concentration at the interface. Counterion concentration 
at the interface {z = D/2) can be, in turn, connected 
with the concentration at the midplane {z = 0), thus 
yielding an alternative form of the osmotic pressure as 



Ii{D) = kBT Cm- 



(19) 
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This latter equation has exactly the same form as in 
the standard PB theory. The only way non-electrostatic 
terms of the free energy influence the osmotic pressure 
is via the solution of eqEl Both forms of the osmotic 
pressure, easlTBland lTQl of course yield exactly the same 
values. 

A typical isotherm Il{D) (in Pascal units) calculated 
using eg s [TBI and [T^ is shown in Figure 2a, and the corre- 
sponding surface charge density in Figure 2b. Note that 
is this example we do not include contributions from hy- 
dration. We discuss the parameter range in the next 
section, and choose the parameters here to be a = —6, 
a = 8 A and x = 12. The isotherm clearly exhibits a 
first-order phase transition from one free energy branch 
at large inter-lamellar separation D to another at smaller 
and with a coexistence region in between. For large 
values of D, D > 64A, most counterions are dissociated 
from surfactants, r^g < 1, and the osmotic pressure fol- 
lows the standard PB theory for (almost) fully dissoci- 
ated surfactants. For D smaller than 39 A, 77^ < 0.1, and 
the isotherm follows another branch, characterized by a 
much smaller surface charge of only about 10% of the 
fully dissociated value. For 39 A < D < 64 A the system 
exhibits a coexistence between two phases. The pressure 
has a plateau and rjs changes from one branch to the 
second. 

The plateau in the osmotic pressure between the two 
stable solutions on Figure 2a is evaluated using the usual 
Maxwell construction. This plateau in 11 (I^) indicates 
a first-order phase transition between the solutions cor- 
responding to two stable branches of the free energy: a 
condensed one with D :^ 39A and a dilute one with 
D :^ 64 A. Although it is hard to fit the phenomeno- 
loRical parameters a, a and x the DDABr results of 
Ref. 7 we believe that the mechanism proposed here 
and the typical results presented in Figure 2 are relevant 
to the experimental system. More details are given in 
Sec. Ill below, where we show how the phase transition 
depends on model parameters. 
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FIG. 2: (a) The osmotic pressure 11 in units of 10^ Pascals (~latm); 
and, (b) the area fraction rjs = a'^a/e of surface charges, as function 
of inter-lamellar spacing D for a = —6, x = 12 and a = 8 A. The 
Maxwell construction gives a coexistence between a phase with D ~ 
39 A and low rjs < 0.1, and another with D ~ 64A and rjs ~ 1. In 
(b) the two coexisting phases are denoted by squares and the dotted- 
dashed line shows the tie-line in the coexistence region. 



C. Added-Salt Case 



After demonstrating the possibility of a phase transi- 
tion for the counterions only, we generalize the results by 
coupling the lamellar system to a reservoir of 1:1 mono- 
valent salt. This generalization is easily achieved within 
our model. The main effect of the salt is to diminish or 
even to eliminate completely the phase transition. 

The main modifications of the model wrought by the 
introduction of salt ions are as follows. First we must con- 
sider ionic concentration profiles independently for posi- 
tive and negative ions: (l)'^{z) and (z). Note that we do 
not distinguish between the negative salt counterions and 
those dissociated from the surfactant. The free energy is 
written, similarly to eq^ in terms of rescaled variables 
with the introduction of two chemical potentials. The 
first, is coupled to the excess amount of co-ions with 



respect to the reservoir, stemming from the integral of ^+ 
in the region between the two plates. The second, 
is coupled with the excess number of counterions, stem- 
ming from the integral of in the region between the 
two plates with the added contribution of surface charges 
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1 — ris- Thus, the total free energy is given by 
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(20) 



Taking now the variation of the above free energy with 
respect to r^g, and we get the Euler-Lagrange equa- 
tions, 



Vs 



exp (/i + a + x^s - ^s) 

exp (-/i" + ^(z)) 
exp - 

3,,-(z)-</>+(z)) 
3 ,^-{z)-^^{z)) . 
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(22) 
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(24) 



The final variation with respect to ys gives the usual elec- 
trostatic boundary condition relating the surface electric 
field with the surface charge density: y^ = {At:Ib / o^)r]s^ 
eq IHl Requiring that the bulk concentration of co-ions 
and counterions matches the reservoir concentration, (/)o, 
where the potential vanishes, = 0, it is easily verified 
from eas l22l24l that 
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(25) 
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sinh {y{z)) = A^^ sinh?/ , (27) 



where 
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D 
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B90 
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(28) 



is the Debye-Hiickel screening length. Inserting eal^into 
1^ leads to a Langmuir-Frumkin-Davis isotherm, now of 
the form: 



Vs 



1 + 0oe"""^'^«+^^ 



(29) 



This equation resembles eq ^21 only that in eq [531 the 

reservoir concentration (j)^ takes the place of (prn. More- 
over, the value of the potential at the midplane is not 
fixed, but rather is determined from (/)o, as are the con- 
centrations, (j)^ = (j)oex.p{'=fym)' 

The PB equation in presence of salt, eql^ depends 
on the electrostatic boundary conditions and has a well- 
known solution expressed via an elliptic integral 4] . The 
first integration of the PB equation, eqEZl from the mid- 
plane position {z = 0) to an arbitrary z gives 



y\z) = V'2 cosh y{z) -2 cosh 
Ad 



(30) 



The boundary condition can be inserted in eqlHOl yielding 



cosh, ys = cosh, yr, 



ttIb 2 



(31) 



A further integration can be written in terms of an elliptic 
function 



z 



dw 



VTcoshi(7^^Tcosh~^ 



and the second boundary condition ys 
be expressed as 



y{z 



D 

Jy^ 



dw 



\J2 cosh w — 2 cosh y^ 



(32) 
L»/2) can 

(33) 



The procedure to solve these equations is similar to the 
one used in the previous section. For given 
X and a, the profiles ^^(2;) and the surface value r]s are 
calculated numerically and inserted into the free energy 
expression, eq EHl Ftot[^=^(^);^s,i^]. Taking the deriva- 
tive of Ftot with respect to or equivalently using the 
contact theorem as for the counterion-only system, the 
osmotic pressure, 11(1)) is obtained. 



n(I)) = /cBT (C++C- -2co) = 



2(/)o). 
(34) 



III. RESULTS 



To better appreciate the role of ionic species in de- 
termining the phase transition we first present results 
for counterions-only and in the absence of hydration in- 
teractions. These results will already clearly show the 
most important feature of the model, namely the possible 
lamellar-lamellar transition. We then predict the effect 
of added salt. Finally, we include the hydration contri- 
bution and compare model results with experiments. 



A. Ion Dependent Lamellar-Lamellar Transition 

Our model contains three parameters: a, x and a. For 
a we use the linear size of the surfactant headgroup, typ- 
ically in the range of 7A to 9A with cross-sectional area 
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FIG. 3: The osmotic pressure isotherm 11(D) for three binding 
strengths: a = —5.85 (dotted-dashes), —5.95 (dashed line), and 
—6.15 (solid line). The other parameters are for non-ideal mixing 
X = 12 and lateral separation a = 8 A. The a = —5.85 and —5.95 
lines show a phase transition, while a = —6.15 does not show one. 
For comparison, the usual PB isotherm with 77^ = 1 (short dashes) is 
also shown. 



a^. For the phenomenological constants a and, in partic- 
ular, for X there is no direct and accurate experimental 
measurement. However, estimates consistent with exper- 
imental data (as discussed in Sees. IIIC and IV), will be 
used hereafter. In our model, positive x presents the pos- 
sibility of lateral phase separation. In addition, for a < 0, 
the counterions tend to stay associated with the charged 
headgroup and reduce the surface charge density. 

In Figure 3 we compare the osmotic pressure isotherm 
Il{D) in the case of no added salt, for three a values (and 
for a constant x = 12 and a = 8 A) with the standard PB 
isotherm (short dashes), for the fully dissociated limiting 
case, rjs = 1. Formally, full dissociation can be achieved 
by setting a ^ 00 in our equations. For the two values of 
a = —5.85, —5.95, the isotherms in Figure 3 show a first- 
order phase-transition in the range of lOA < D < 50 A. 
The phase transition is from a dilute and highly charged 
Lq, lamellar phase (large D and r]s ^ 1), to another Lq,/ 
phase that is more condensed and less charged (small D 
and 1). As a increases, this phase transition shifts 

to higher values of H and lower values of D. 

For large r]s is very close to one, and the osmotic 
pressure isotherm closely follows the PB result; for large 
D the solution of our model remains essentially on the 
PB branch of the free energy, characterized by 77^ :^ 1. 
For small values of inter-lamellar spacing, the values 
of osmotic pressure for all three values of a are again very 
similar, with small and slowly varying 77^, as in Figure 2b. 
Here, the system essentially remains on the associated 
branch of the free energy characterized by 7^5 <C 1. 

The lowest a isotherm (solid line, a = —6.15) shows no 



FIG. 4: The osmotic pressure isotherm 11(D) for three values of 
the surfactant headgroup separation: a = 7.5 A (short-dashes), 8 A 
(dashed line), and 8.5 A (solid line). The other parameters are a = —6 
and X = 12. A phase transition is seen for a = 8 A and 8.5 A, but not 
for a = 7.5 A. 



transition. The counterions are almost fully associated 
in this case for the entire range of leading to a lower 
value of pressure for all I^'s. Note that all three osmotic 
pressure isotherms as well as the PB one have the same 
limiting behavior for D ^ 00. This is quite accurately 
described by the Langmuir form of the osmotic pressure, 
as applied to the counterion-only case 3 : 



n(D) 



ksTn 1 
21b 



(35) 



which does not depend on the value of the surface change 
(or, equivalently on r]s). 

In Figures 4 and 5 we show the effects of the variation 
of linear size a and lateral interaction x^ respectively. 
We have chosen the parameter range to show isotherms 
without a phase transition (small a or x) as well as those 
showing the transition (higher values of a or x) each 
of the figures. The main features are the same as in Fig- 
ure 3. The large D region represents highly dissociated 
lamellae (strongly charged), while the phase transition 
(when it exists) can be seen for small D at higher lamel- 
lar density. Increasing a moves the transition point to- 
wards higher values of osmotic pressure or equivalently 
lower values of D until it eventually disappears. Increas- 
ing X has the same effect. Note also that the pressure at 
low D is diminished on increase of a. This effect can be 
understood by recalling the relation a = erjs/d^i so that 
for the same 7^5, larger a corresponds to smaller surface 
charge density a. 

Our results are summarized in Figure 6 where the (a, 
x) parameter space is separated by a solid descending 
line Xc(<^) into two regions (for fixed a). The upper re- 
gion (large Xi large a) corresponds to isotherms with a 
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D[k] 

FIG. 5: The osmotic pressure isotherm 11(D) for three x values: 
X = 11 (short dashes), 12 (dashed line), and 12.25 (solid line). The 
other parameters are a = —6 and a = 8 A. The phase transition is 
seen for the two larger values of x- 

phase transition (and is designated as "two phases" on 
the figure). Below that region (smah smah a) the 
isotherms show no phase transition (designated as "one 
phase" on the figure). The degree of counterion disso- 
ciation varies in this region from very smah values to 
values 775 < 0.8 for finite values of a and x- The PB re- 
sult of 775 = 1 is reached only asymptotically as a ^ 00. 
The line represents the continuous line of critical point 
in the (x,Cif) plane. The region between the solid and 
dot-dashed lines corresponds to jumps in D at the tran- 
sition of more than 3 A. Above the dot-dashed line, the 
behavior at I) > 3 A is described by the usual PB solu- 
tion because the transition occurs at unphysically small 
D values. 

The inset to Figure 6 corresponds to variations of x 
and a shown in Figure 4 and 5, respectively. As x oi" <^ 
are lowered, the single-phase region is eventually reached. 
Figure 6 is plotted for a = 8 A. As a decreases, the two- 
phase region shrinks and the one-phase region expands. 

The dashed line in the figure is an analytic calculation 
which gives the following approximate form of x*(<^) 

X*c^-2a + 21n(^^^ -4. (36) 

To derive this result we assume that the phase transition 
occurs at large D. Using the asymptotic large D behav- 
ior, we compare the free energy of 775 :^ 1 with 775 <c; 1 
and determine the transition point as a function of x for 
given a and a. As can be seen by comparing the analytic 
(dashed) line with the full numerical solution (solid line), 
the approximation is good for small a < 0. For a > — 3 
the assumption of a transition at large D breaks down 
and the approximated x* starts to deviate considerably 
from the numerically calculated Xc(<^)- 
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FIG. 6: Phase diagram in the {a, x) plane for a = 8 A. In the region 
below the solid line Xc((^), the system does not exhibit a phase tran- 
sition ("one phase"). The dashed line is the analytical approximation 
of X* — — 2q; + 2 ln(27r/ b /a) — 4. The inset is a blow-up of the region 
around a ~ —6 and x ~ 12, showing the parameters used in figures 
3 and 5. Square symbols correspond to isotherms in Figure 3, while 
diamonds correspond to those in Figure 5. The region lying between 
the full and dot-dashed lines in the two-phase region corresponds to 
transitions in D that are larger than ?>h. 

B. Added Salt: Vanishing of the Transition 

The effect of added salt was treated in Sec. lie. 
The salt is characterized by the Debye-Hiickel screening 
length Ad, and screens electrostatic interactions. As the 
amount of added salt increases, \d decreases, and the 
phase transition observed in the absence of salt becomes 
gradually less pronounced until it is finally wiped out 
completely 45]. This is clearly seen in Figure 7. In 7a, 
three osmotic pressure isotherms are shown. A plateau 
(first-order phase transition) is seen for the two lower 
amounts of salt, C5 = 10 and 30 mM, while the transition 
disappears for higher amounts of salt, C5 = 50 mM. 

One can also see how the phase transition is first 
shifted towards the high D low H values, and then (for 
C5 ~ 40 mM) completely disappears. The overall de- 
crease in n(I)) as the amount of salt increases, is due 
to the increased screening, and is present also in the sim- 
ple PB theory. In 7b the jump in rjg is shown for C5 = 10 
and 30 mM, while it vanishes for higher amounts of salt, 
C5 = 50 mM, in accord with the isotherm behavior. 

It is instructive to follow the change of the transition 
pressure Iltr as salt is added to the system, which is re- 
lated to the difference in volume AV and number of ions 
AA^ in the two phases in a Clausius-Clapeyron-like equa- 
tion: 

dcb Cfc AT/' ^ ' 

Remarkably, we find an almost linear dependence of Iltr 
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FIG. 7: Effect of added salt. In (a) U{D), in (b) r]s{D) are plot- 
ted as a function of added salt concentration c^. In (a) and (b): 
= 10 mM (solid line), 30 mM (dashed line), 50 mM (short dashes). 
In (c) salt concentration varies continuously over the range of ex- 
istence of the transition, < < 40 mM. Other parameters are 
ck: = —5, X = 10.19, and a = 8 A. As increases, screening becomes 
more important. The entire osmotic pressure isotherm 11(D) decreases 
in magnitude and the phase transition region diminishes and shifts 
towards higher D values. Note that for highest salt concentration, 
Cb = 50 mM, the phase transition has vanished. In (b) the coexisting 
values of the two phases are denoted by a square (c^ = 10 mM) and by 
a circle (c^ = 30 mM). The dotted-dashed lines are the corresponding 
tie-lines. 



in the whole range of C5, starting with the transition 
pressure at no added salt and leading eventually to the 
loss of transition for sufficiently high salt concentrations, 
C5 ^ 40 mM, see figure 8. In eqEI dUtr/dci) corresponds 
to an added work of 2 ksT due to the exclusion of ions 
acting on the volume change at the transition. 

It is also interesting to follow the change in Iltr with 
temperature. Experimentally, both an increase of Iltr 
with T (at lower T) and a decrease (at higher T) have 
been observed In contrast, if we assume that x 
a are independent of T in the model, we find that Iltr 
monotonically, almost linearly, increases with T. This 
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FIG. 8: Effect of added salt on the transition pressure. Iltr{cb) is 
plotted as a function of added salt concentration c^. The dependence 
of the transition pressure on the salt concentration in the region where 
the transition exists is linear, with a slope of 7.9 x 10~^^ J. 



suggests that in a more refined model, the parameters x 
and a should be taken as functions of the temperature 
rather than simple constants. For example, if we assume 
that a{T) and x{T) vary as 1/T, while the ratio a/x is 
kept constant, we find that Iltr decreases monotonically 
with T. 



C. Relating to DDABr/DDACl Osmotic Pressure 
Experiments 

The experimentally observed difference between 
DDABr and DDACl for the different halides can easily 
be rationalized within our model by different values of a 
and/or for the different ions. This is indeed reason- 
able since experiments show that larger halide ions have 
an added affinity even for neutral lipids [23, • 

The tendency of ions to preferentially partition into the 
hydrocarbon-water interface is most often reported in 
terms of an effective binding interaction that acts in addi- 
tion to the repulsive electrostatic force, expected for ions 
interacting with low dielectric material. These differences 
in binding affinity would translate into a different value 
of a within our model. In these terms, experiments show 
that |a| is larger for bromide by one to four ksT more 
than for chloride, and iodide is at least an order of mag- 
nitude larger than those mmiH^]. 

With our model assumptio ns, we can now try to fit 
the experimental data in Ref. [l^ using the same (small) 
amount of added salt as in the experiment, i.e. C5 = 
0.5 mM. The fit to the DDABr and DDACl lamehar sys- 
tems are shown in Figure 9. We will first address the 
fit to the simpler case of DDACl that does not show in 
experiments a liquid-liquid coexistence, and then discuss 
DDABr, where the liquid-liquid coexistence is clearly dis- 
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FIG. 9: Fit to the experimental osmotic pressure isotherm n(D) of 
Ref. 3 log-log scale. The diamonds and squares are the data 

points for DDABr and DDACI, respectively, reproduced from Ref. IT4t . 
In (a) the solid line is the best fit of the model to the phase transi- 
tion seen for DDABr with a = -7.4, x = 14.75 and a = 8 A. The 
fit also includes a hydration contribution (parameters for the form: 
nhyd = no exp(-D/Ahyd), with typical values: Uq = 2.37 • 10^ Pa 
and Ahyd = 1.51 A). This contribution is particularly important at the 
low D region of the DDABr isotherm. A small amount of salt is added 
in the fit as in the experiment (c^ = 0.5 mM). The dot-dashed line 
is the fit to the DDACI (no transition). All parameters here are the 
same, except a = —3.4 (one-phase region of Figure 6). In (b) the 
same model fits, showing only the electrostatic contributions to the 
force, Fv + Fg, that are responsible for the apparent phase transition. 
The fit to both data sets is good. However, the steep increase in 
pressure of the DDABr seen in the low D range (D < 10 A) cannot 
be reproduced with these forces alone. 



cerned. The DDACI data was fitted in Figure 9a us- 
ing a = -3.4, X = 14.75 and a = 8A. The DDACI 
data points, represented by squares, are reproduced from 
Ref. The value chosen for x is higher, yet close to 

Xc(<^) (See Figure 6). 

In comparing to experiments, it is important to realize 



that for small spacings, interactions of non-electrostatic 
origin tend to dominate the osmotic equilibrium 0, 0, 
0, 0, Even for highly charged systems, hydra- 

tion interactions acting at very short range, practically 
independent of the charging equilibrium at the surface, 
invariably dominate at separations of D < 10 A are cer- 
tainly among the most important to consider. There- 
fore, we add the hydration interaction Fhyd as an ide- 
alized separable term to the overall interaction energy: 
Ftot + ^hyd = Fy -\- Fs -\- Fhyd • couformity with hydra- 
tion interactions between lipids and between other 
macromolecules we use the phenomenological form 
of an exponential interaction with a salt-independent de- 
cay length, Ahyd [Hoj 

^hyd = no Ahyd exp(-i:>/Ahyd). (38) 

To account for the hydration interaction, we fit with IIo = 
2.37 X 10^ Pa and Ahyd = 1-51 A. These values are similar 
to those found from fits to experimental inter-lamellar 
spacings of lipids 0, 0, 'si] . 

Note that our extended PB model predicts a transi- 
tion at these values of a and x- However, this phase 
transition occurs at very small inter-lamellar separations 
of D < 3A, indicating that only at such low D values 
are the ions associated (small r]s). For higher D values, 
the ions dissociate and the model follows the regular PB 
osmotic pressure. Values of I) < 3A are well below the 
validity limit of the model; for such low I^'s the computed 
osmotic pressure isotherm can not be considered realis- 
tic. For higher the surface charge density, 775, is almost 
one and the system exhibits the usual Poisson-Boltzmann 
isotherm, with screened electrostatic interactions at large 
D. We remark that the fit in the high D region is very 
good, and surprisingly, also persists to small D (although 
there are very few data points below I) = 30 — 40 A) . 

The DDABr was fitted in Figure 9a using a different 
a = —7.4 but the same x = 14.75 and a = 8 A, as well 
as the same hydration interaction. The value of x was 
picked near Xc so that the transition would occur at a 
value of D in the range of 10 to 40 A. In this range our 
model is still valid and nicely reproduces a coexistence 
between the two liquid lamellar phases La and Lc^/. 

To better appreciate the role of hydration, figure 9b 
shows the model fits that exclude the hydration contribu- 
tion to the free energy. At smaller separations, D < lOA, 
the fit deviates considerably from the experimental os- 
motic pressure isotherm. While electrostatic forces can 
account for the transition point itself, they can definitely 
not by themselves reproduce the sharp rise in the osmotic 
pressure in the undissociated branch following the phase 
transition, compare Figure 9a. 

Juxtaposing Figures 9a and 9b, it is evident that 
adding the hydration contribution affects mainly the 
undissociated branch at small inter-surface separations, 
since at large spacings electrostatic PB contributions 
completely overwhelm the much shorter ranged hydra- 
tion. The hydration interaction mainly affects the small 
D region of both DDACI and DDABr. In fact, while 
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the magnitude of hydration forces may not differ greatly 
from other (electrostatic) acting forces, the rate of force 
decay (slope) at such low spacings indicates that hydra- 
tion is always acting. One can observe this in the context 
of lipids 48] as weh as DNA 52]. Fits to both DDACl 
and DDABr should, therefore, contain hydration contri- 
butions to the free energy. 

Importantly, as can be realized from Figure 6, there 
are different ways to cross the transition line, Xc{<^)^ to 
witness the liquid-liquid coexistence. Another possible 
choice of parameters could thus assume the variation in 
X instead of a. For example, the data can be equally 
well fit if we fix a = —7.4 and fit the DDABr data with 
X = 14.75, and the DDACl with x = 18- 

Finally, we address the iodide analog, DDAI, that in 
experiments appears not to form a swollen liquid phase, 
remaining in the condensed liquid lamellar phase for the 
whole range of osmotic pressure values 0. In the lan- 
guage of our model DDAI has either large \a\ or x val- 
ues; it never enters the PB branch but instead remains 
in the undissociated branch for all spacings D. This cor- 
responds to the "one phase" region of figure 6. Indeed, 
\a\ for iodide could be estimated as yet larger than in 
the case of DDABr, probably by up to an order of mag- 
nitude |2^ [2^ . The system is confined to remain on the 
undissociated branch for this counterion: the repulsion 
due to electrostatic interactions and entropy of counter- 
ions is not strong enough to combat the (van der Waals) 
attractions. These attractions that are strong enough 
to hold the system in the secondary free energy mini- 
mum, effectively prevent the swelling observed for the 
other two counter ions. We could thus establish a ranking 
in the model's parameter space, where DDACl, DDABr 
and DDAI would make a monotonic Hofmeister-like se- 
ries either in the values of \a\ or x- The 3rd parameter, 
a, is not changed in the fit because it is taken as the size 
of the same DDA+ headgroup. 



IV. DISCUSSION 



The model presented in this paper combines electro- 
static and non-electrostatic interactions between charged 
surfactant bilayers. The non-electrostatic part of the in- 
teractions is accounted for using two phenomenological 
parameters characterizing the strength of the counterion- 
amphiphile interaction, a, and lateral amphiphile- 
amphiphile interactions, x^ the charged dissociable 
surface. By choosing a to be negative, corresponding to 
favorable adsorption energy, and x to be positive, cor- 
responding to a net attractive interactions between like- 
species, promoting lateral segregation between the dis- 
sociated and non-dissociated surfactants on the surface, 
we are able to qualitatively explain the experimentally 
observed Lq, Lq,/ lamellar-lamellar phase transition. 



A. Origin of Phase Transition 

An abundance of experimental verification [25I . [25 [2^ 

28] indicates that the different halides interact differently 
with lipids, forcing us to recognize the existence of non- 
electrostatic interactions at work. In our model, this 
preferential interaction is represented by a. By reducing 
the effective layer charge density, a favorable preferen- 
tial interaction of ions to the interface acts to lower the 
pressure at any given spacing D. However, a alone can- 
not account for the abrupt jump in D associated with 
a first-order phase transition. In fact, in the absence of 
X, pressure isotherms for any value of a are smooth (no 
phase transition). 

Our model offers a natural extension of the PB theory 
with the Langmuir-Frumkin-Davies adsorption theory as 
applied to simple counter ions. The large difference in be- 
havior of the halide ions is modeled by the parameters 
a, X cind a. The interaction parameters a and x ^^c- 
essarily involve contributions from changes in hydration, 
solvation and desolvation, of lipid headgroups and their 
counterions. The model, therefore, underscores the im- 
portant role of water structuring around ions devolved in 
the bulk versus at the interface. 

The salient feature of our model is the first-order tran- 
sition in the osmotic pressure isotherm from an almost 
completely dissociated state (highly charged and swollen, 
PB branch) at large interlayer separations, to an almost 
neutral, weakly dissociated, state approaching bilayer 
contact (condensed, undissociated branch). The exis- 
tence of this transition depends on the values of both 
phenomenological parameters, but it is present over an 
extended region of phase space. The PB branch of the 
osmotic pressure isotherm is not much different from the 
standard PB theory with complete dissociation, both 
with or without added salt. On the other hand, the 
undissociated branch is characterized by a large atten- 
uation in the magnitude of the osmotic pressure for a 
given inter-lamellar spacing, being about two orders of 
magnitude smaller then in the PB case. 

In our model the Lq, L^' transition in the inter- 
lamellar spacing is coupled to a lateral first-order phase 
transition of the r]s order parameter. This is a direct con- 
sequence of the coupling between inter-lamellar electro- 
static degrees of freedom (mean electrostatic potential, 
mean ion density) and the surface non-electrostatic de- 
grees of freedom as quantified by the phenomenological 
parameters a and x- The ensuing liquid-liquid lamel- 
lar phase transition is thus not only from one state of 
the lamellae where a larger fraction of the amphiphiles 
are charged to another state where they are less charged, 
but also from a state where the inter-lamellar forces are 
by and large electrostatic in nature, to a state where they 
are dominated by hydration. While in many experimen- 
tal systems this transition is smooth and gradual, it is 
quite pronounced and discontinuous in the system stud- 
ied here and in Ref. |T^ . 
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B. Relating Model Parameters and Molecular 
Interactions 

We propose that non-ideal mixing between counterion- 
associated and dissociated surfactants can be responsible 
for an in-plane transition, which, in turn, is coupled to 
the bulk transition. This proposed non-ideality is repre- 
sented in our model by as is sometimes used to report 
on lipids showing phase transitions following changes in 
pH Is^- While at present direct experimental verifica- 
tion and estimates for the proper x values are lacking, we 
propose that conformational changes of lipid headgroups, 
ions and water structuring induced by the adsorbing ion, 
together with an added van der Waals interaction be- 
tween neutralized surfactant complexes can lead to sig- 
nificant demixing. Furthermore, because larger ions are 
expected to perturb the lipid- water interface to a greater 
extent, it is reasonable to expect that the value of x will 
follow a similar ranking to the binding of ions to the in- 
terface, represented by a. We note, that the x values 
needed to observe a phase transition, typically ^ lO/c^T, 
are quite high. These high values are needed to overcome 
the electrostatic repulsion between like-charged lipids in 
this unscreened, highly charged system. The source of 
this lipid demixing energy (our x parameter) could be 
associated with mismatch of headgroup-headgroup in- 
teractions, such as hydrogen bonding between neutral 
lipids, water-structuring forces, or non-electrostatic ion- 
mediated interactions between lipids across two apposed 
bilayers for small inter-lamellar separations. 

The parameter models the area per headgroup on 
the membrane plane. It is a function of several molecular 
interactions and, in principle, can be determined varia- 
tionally. In Ref. 143 area/headgroup was found to 
vary in a non-trivial fashion, from a larger value in the 
condensed lamellae to a smaller one in the dilute lamel- 
lae. The forces determining the area per surfactant are 
as yet unknown. Therefore, in the model we have not al- 
lowed for changes in area per surfactant, but note that it 
is not inconsistent to assume that the area per headgroup 
differs for the neutral vs. charged surfactant. The expan- 
sion of lipid area upon condensation, contrary to what is 
typically observed in phase transitions of lipids, could be 
evidence for direct attraction between Br~ ions and the 
lipid hydrocarbon core, as suggested previously by Nin- 
ham et al 0. This point deserves further investigation. 
Interestingly, however, as seen in Figure 4, the model 
predicts that lipids with larger area per headgroup can 
show a larger transition gap, as was found experimentally 
in the case of GM-3 ganglioside with weakly adsorbing 
counterions j5^ . 

Another point that merits further investigation is the 
dependence of the isotherms on temperature. Both a{T) 
and x(^) are complex functions of the temperature, with 
specific dependence that cannot be obtained from our 
model. A change in temperature affects the values of 
a{T) and x{T) and can change the plateau pressure 
values, as was measured and reported in Figure 9 of 



Ref. A better understanding of a{T) and x(^) ^ay 
offer an explanation to the non-monotonic behavior of 
the plateau pressure as function of temperature 14]. 

Finally we point out that the lack of direct experi- 
mental evidence at this time, particularly for x? hmits 
out predictions to be mainly qualitative. More specifi- 
cally, we can offer only a qualitative explanation for the 
strong difference in the behavior for different counterions, 
namely: no transition, transition and no stable swollen 
lamellar phase for the DDACl, DDABr and DDAI am- 
phiphilic systems, respectively. Therefore, the fit to the 
data points shown on Figure 8 should be regarded as a 
tentative explanation of the mechanism behind the ob- 
served phase transition. For example, we cannot estab- 
lish whether the main difference in ionic interactions with 
the surface are properly characterized by the value of a 
as opposed to x- 



V. CONCLUDING REMARKS 

Previously, phase transitions in DDABr lamellar sys- 
tems have been theoretically attributed to either an 
ion-dependent van der Waals attraction between layers 
6, 54], or to a strong-coupling effect between adsorbed 
ions, expected for surfaces with high charge density [HEj . 
Here, we have shown that it is possible to account for the 
phase transition assuming a non-electrostatic interaction 
between ion-dissociated and ion-bound surfactants. We 
suggest that this interaction is ion specific and, hence, 
we offer an explanation for the different behavior seen 
for the three halide counterions. 

The large phenomenological parameters we have found 
in our own fits of the data (see Fig. 9) as well as the 
large energetic terms assumed in the other approaches 
^, 54, 55], all indicate that substantial attraction nec- 
essarily acts to overcome the electrostatic repulsion be- 
tween surfactants. The molecular origin of this large en- 
ergy has not yet been determined, and further experi- 
mental verification of the different phenomenological pa- 
rameters is required. It will also be of interest to see if 
the lateral phase transition underlying the swelling tran- 
sition and its dynamical evolution from one lamellar state 
to the other, can be observed directly in experiments. 

We hope that in follow-up studies, a more microscopic 
approach will be able to shed light on the origin of the 
phase transition in these charged lamellar systems and 
how they relate to specific molecular details. 
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